RADII MINIMAL PROJECTIONS OF POLYTOPES AND 
CONSTRAINED OPTIMIZATION OF SYMMETRIC POLYNOMIALS 

RENE BRANDENBERG AND THORSTEN THEOBALD 

Abstract. We provide a characterization of the radii minimal projections of polytopes 
onto j-dimensional subspaces in Euclidean space E™. Applied on simplices this char- 
acterization allows to reduce the computation of an outer radius to a computation in 
the circumscribing case or to the computation of an outer radius of a lower-dimensional 
simplex. In the second part of the paper, we use this characterization to determine the 
sequence of outer (n — l)-radii of regular simplices (which are the radii of smallest en- 
closing cylinders). This settles a question which arose from the incidence that a paper by 
Weifibach (1983) on this determination was erroneous. In the proof, we first reduce the 
problem to a constrained optimization problem of symmetric polynomials and then to an 
optimization problem in a fixed number of variables with additional integer constraints. 



I. Introduction 

Let Cj tn be the set of all j-dimensional linear subspaces (hereafter j-spaces) in n- 
dimensional Euclidean space E n . The outer j -radius Rj(C) of a convex body C C E" is 
the radius of the smallest enclosing j-ball (j-dimensional ball) in an optimal orthogonal 
projection of C onto a j-space J G £j, n , where the optimization is performed over L^ n . 
The optimal projections are called Rj -minimal projections. In this paper we show the 
following results: 

Theorem 1. Let 1 < j < n < m and P = conv{u^, . . . ,v^} C E n be an n-polytope. 
Then one of the following is true. 

a) In every Rj-minimal projection of P there exist n + 1 affinely independent vertices of 
P which are projected onto the minimal enclosing j -sphere. 

b) j > 2 and Rj(P) = Rj-\{JP H H) for some hyperplane H = aS{v^ : i E 1} with 
I C {l,...,m}. 

If j = 1 or if P is a regular simplex then always case a) holds. Moreover, the number v 
of affinely independent vertices projected onto the minimal enclosing j -sphere is at least 
n — j + 2 and there exists a (v — 1 )-flat F such that Rj{P) = i2j+^_„_i(P D F). 
The bound n — j + 2 is best possible. 
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Theorem ^ allows to reduce the computation of an outer radius of a simplex to the 
computation in the circumscribing case or to the computation of an outer radius of a 
facet of the simplex. 

Using this theorem, the second part of the paper shows the following result on the outer 
(n — l)-radius, which is the radius of a smallest enclosing cylinder. 



Theorem 2. Let n > 2 and T™ be a regular simplex in E n with edge length 1. Then 

R n -i{Ti) = 



n-1 



2(n+l) 
2n-l 



2y/2n(n+l) 



if n is odd, 
if n is even. 



The case n odd has already been settled independently by Pukhov [TQ and WeiBbach [22] , 
who both left the even case open in their papers. Pukhov's results also determine Rj(T™) 
for j < n. There also exists a later paper on i? n _i(T n ) for even n j2S], but as pointed out in 
(Sj the proof contained a crucial error. 1 Thus Theorem |21 (re-) completes the determination 
of the sequence of outer j-radii of regular simplices (see also |2j). 

Studying radii of polytopes is a fundamental topic in convex geometry. Motivated 
by applications in computer vision, robotics, computational biology, functional analysis, 
and statistics (see jH] and the references therein) there has been much interest from the 
computational point of view. See [3J HI] for exact algebraic algorithms, P [2D 121] for 
approximation algorithms, and [U |H] for the computational complexity. Reductions of 
smallest enclosing cylinders to circumscribing cylinders are used in exact algorithms as 
well as for complexity proofs (see, e.g., (21 Theorem 1] and [SJ Theorems 5.3-5.5]), and 
have previously been given only for j £ {1,^} as well as for dimension 3. Theorem^ 
generalizes and unifies these results. 

Here, we use Theorem ^ to reduce the computation of the outer (n — l)-radius of a 
regular simplex to the following optimization problem of symmetric polynomials in n 
variables: 



s.t. 



n+l 




min ^2 sf 




i=i 




n+l 




£*? 


= o, 






n+l 






= i, 


i=i 




n+l 




E s i 


= 0. 


i=i 





Based on exploiting the symmetries, we solve (jl.lj) for any n by reducing it to an opti- 
mization problem in six variables with additional integer constraints. 



After 3 had been completed, Bernulf WeiBbach suggested to work jointly on a new proof for the even 
case. Unfortunately, he died on 9th June 2003. 
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The paper is structured as follows. In Sectional we introduce the necessary notation. 
Section |3] gives the proof of Theorem ^ Section H] contains the derivation of the opti- 
mization problem and the proof of Theorem El In Section |3] we analyze the different 
difficulty of the even and the odd case of (jl.lj) from the viewpoint of the Positivstellensatz. 

2. Preliminaries 

Throughout the paper we work in Euclidean space E n , i.e., ¥L n with the usual scalar 
product x ■ y = Y^i=i x iVi an d norm = [x ■ x) 1 / 2 . B n and S n_1 denote the (closed) 
unit ball and unit sphere, respectively. For a set A C E", the linear hull of A is denoted 
by lin(y4), the affine hull by aS(A), and the convex hull by conv(A). 

A set C C E n is called a body if it is compact, convex and contains an interior point. Ac- 
cordingly, we always assume that a polytope ? C E" is full-dimensional (unless otherwise 
stated). Let 1 < j < n. A j-flat F (an affine subspace of dimension j) is perpendicular 
to a hyperplane H with normal vector h if h and F are parallel. For p,p' e E n and 
subspaces E G Cj jU , E' e Cf, n , a j-flat F = p + E and a j -flat F' = p' + E' are parallel 
if EUE' = \in{EUE'). 

A j- cylinder is a set of the form J + pW 1 with an (n — j)-flat J and p > 0. For a 
body C C E n , the outer j-radius Rj{C) of C (as defined in the introduction) is also the 
radius p of a smallest enclosing j-cylinder of C. It follows from a standard compactness 
argument that this minimal radius is attained (see, e.g., |2J). Let 1 < j < k < n. If 
C C E" is a compact, convex set whose affine hull F is a fc-flat then Rj(C) denotes the 
radius of a smallest enclosing j-cylinder C relative to F, i.e., C = J' + Rj(C')(M n PI F) 
with a (k - j)-flat J' C F. 

A simplex S := conv{t> . . . , t>( n+1 )} (with v^ 1 ', . . . , t>( n+1 ) e E n affinely independent) 
is regular if all its vertices are equidistant. Whenever a statement is invariant under 
orthogonal transformations and translations we denote by T n the regular simplex in E n 
with edge length \/2. The reason for the choice of a/2 stems from the following embedding 
of T n into E n+1 . Let H n a = {x E E n+1 : J27=i x i = «}■ Then the standard embedding T n 
of T n is defined by 

T n := conv {e (i) G E n+1 : l<i<n + l}c^™, 

where eW denotes the z-th unit vector in E n+1 . By l S n_1 := S n fl 7ig w e denote the set of 
unit vectors parallel to 7i™. 

A j-cylinder C containing some simplex S is called a circumscribing j-cylinder of S if 
all the vertices of S are contained in the boundary of C. 

3. Minimal and circumscribing j-cylinders 

It is well known that the (unique) minimal enclosing ball B (i.e., the minimal enclosing 
n-cylinder) of a polytope P C E n may contain only few vertices of P on its boundary 
(e.g., two diametral vertices) PU p. 54]. However, in cases where less than n + 1 vertices 
of P are contained in the boundary of B, it is easy to see that there exists a hyperplane 
H such that P fl bd(B) C H and the center of B is contained in H. Then the smallest 
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enclosing ball of P and the smallest enclosing ball of P D H relative to H have the same 
radius. 

In Theorem 1.9] the following characterization for the minimal enclosing 1-cylinder 
(two parallel hyperplanes defining the width of the polytope) is given: 

Proposition 3. Any minimal enclosing 1-cylinder of a polytope PcE™ contains at least 
n + l affinely independent vertices of P on its boundary. 

We give a characterization of the possible configurations of minimal enclosing j-cylinders 
of polytopes for arbitrary j, unifying and generalizing the above statements. 

Lemma 4. Let P = conv{t> . . . , v ^ m '} be a polytope in E n ; 1 < j < n — 1, and J be an 
(n — j)-flat such that C = J + Rj(P)M n is a minimal enclosing j-cylinder of P. Then for 
every I C {1, . . . ,m} such that {% : i>W £ bd(C)} C / and Hi := &S{v^' : i G 1} is of 
affine dimension n — 1, J is parallel to Hi. 

Proof. Suppose that there exists a hyperplane H := Hi of this type with J not parallel 
to H. Let n := \{v^ E H : 1 < i < m}\. Without loss of generality we can assume 
H = {x E E n : x n = 0} and I = {v^\. . .,v^}. Hence, v& +1 \ . . . , w (m) 0#Ubd(C). 

First consider the case that J is not perpendicular to H. Let p, s^, . . . , s^ - ^ e E n 
such that J = p + lin{s^, . . . , s^ n ~^}. Since, by assumption, J is not parallel to H, 
we can assume p = G J C\ H, Sn^ = ■ ■ ■ = s„ J ' — and si" ^ > 0. For every 
< G (0, s { n~ 3) ) and s' := (s^, • • • , sfcf, <) G E n let J' = p + \in{s^, . . . , s^" 1 ), s'}. 
Geometrically, J' results from J by rotating J towards the hyperplane H in such a way 
that the orthogonal projection of J onto remains invariant (see Figure [1J. 




Figure 1. For n = 3 and j = 2 the figure shows how the underlying 
flat J of the j-cylinder C is rotated towards its orthogonal projection onto 
the plane H. The distances between the vertices 1 < i < n, and the 
j-cylinder axis are not increased, and decreased if G" K. 
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Since J and H are not perpendicular we obtain J ^ J', and because . . . , v ^ G H 
that 

(3.1) dist(v w , J') < dist(i; {i) , J) , 1 < i < n , 

where dist(-, •) denotes the Euclidean distance. In (|3.1jh "<" holds whenever i;W ^ K : = 
J^DH. Obviously, dim(F) = j — 1. If none of the lies in nbd(C) then, by choosing 
s' n sufficiently close to sl™ - ^, all vertices of P lie in the interior of C = J' + Rj(P)M n , a 
contradiction to the minimality of C. Hence, there must be some vertex of P in Kf]bd{C). 

Let k := \{v^ G K(lbd(C) : 1 < i < m}|. By renumbering the vertices we can assume 
that v {1 \ . . .,v(~ k) G K n bd(C). Let F := convjV 1 ), . . . , v {Jt) ) and k := dimF. 

Now assume Ffl J = 0. We have shown above that for sufficiently small s' n the rotation 
from J to J' keeps all vertices within the j-cylinder C and t)W,,..,w' fc ' are the only 
vertices on bd(C'). Let J" be a translate of J' with dist(J",F) < dist(J', F), and J" 
sufficiently close to J' to keep . . . ,t>( m ) within the interior of C" = J" + Rj(P)M n . 

Then all vertices of P lie in the interior of C", again a contradiction. 

It follows that F n J ^ 0, and since F C K = J ± n H that F fl J = {p} = {0}. 
Since dist(p, v ^) = Rj(P) for all i & {1, . . ., fc} and since p G F, it follows that p is the 
unique center of the smallest enclosing fc-ball of F. Now let J'" result from J' by rotating 
J' around the origin towards a direction in IR n \ (Ui=i(' y ^)~ L ) > This set of directions is 
nonempty, since every subspace in the union is only of dimension n — 1. For i G {1, . . . , k} 
the property dist(uW, J) = dist(uW, J') = dist(t)W,p) implies dist(t> W, J'") < dist(i>W, J'). 
Again, by keeping the rotation sufficiently small, t>( fc+1 ), . . . , t>( m ) remain in the interior of 
C" = J'" + Rj(P)M n . Now, all vertices lie in the interior of C", once more a contradiction. 

Finally, consider the case that J is perpendicular to H. Then J fl if is an optimal 
(n — 1 — j)-flat for the j-radius of P fl H (taken in (n — l)-dimensional space). However, 
it is easy to see that in this case any small perturbation J' of J with J' fl H = J fl H 
keeps within the j-cylinder, not increasing the distances of all the other 

vertices to the new (n — j)-flat. Indeed, the case J perpendicular to H describes a local 
maximum. So the same argument as in the non-perpendicular case applies to show a 
contradiction. □ 

Corollary 5. Let S = conv{?;W, . . . ; ^( n + 1 )} he a simplex in E n and be the facet of S 
with «W g Letl < j < n-l and J be an (n-j)-flat such that C = J + R j (S)M n is a 
minimal enclosing j- cylinder of S. Then J is parallel to every for which, v® G'bd(C). 

Lemma 6. Let P = convji/ 1 ), . . . ,-i/ m )} be a polytope in E n ; 1 < j < n, and J be an 
(n — j)-flat such that C = J + Rj(P)M n is a minimal enclosing j-cylinder of P. If there 
exists a hyperplane Hj = aS{v^' : i G 1} which is parallel to J, then one of the following 
holds: 

a) There exists a vertex i)W ^ Hj that lies on the boundary of C; or 

b) j > 2, J c Hi, and Rj(P) = Rj^P n Hi). 

Proof. By Proposition El for j = 1 always a) holds; so let j > 2, and suppose neither a) 
nor b) holds. Since b) does not hold there exist (n — j) -fiats parallel to J and closer to 
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nice 



Hi, and since a) does not hold, for any such (n — j)-flat J', such that all vertices £ Hj 
stay within C, the distances from the vertices i e I, to J' are strictly smaller than 
their distances to J. Hence C cannot be a minimal enclosing cylinder. □ 

In the case that P is a simplex, the proof can be carried out more explicitly: 
Let p( n+1 ) be the facet of P not including the vertex t>( n+1 ). Suppose that J is parallel 
to P(" +1 ), that P( n+1 ) C H := {x G E n : x n = 0}, and that > 0. Let p e J. S 

^ > it follows p n > and obviously 

(3.2) R 3 {P)>v^- Vn . 
On the other hand, since J is parallel to p( n+1 \ 

(3.3) Rj(PT = RU( p(n+1) )+ Pi- 
Let 

(v n n+1) ) 2 - ( P («+i)) 



(3.4) p* 



be the unique minimal solution for p n to (|3.2j) and ()3.3|) . Due to p n > 0, we obtain 
p n = max{0,p*}. Now, we see that case a) holds if p n = p* n and case b) if p n = 0. 



Statements almost complete the proof of Theorem ^ If the number v of affinely 
independent vertices of P lying on the boundary of C is at most n, it follows from Lemma 
E] and El that case b) of Theorem ^ must hold. Moreover, if v < n — 1 we can again apply 
these lemmas on the lower-dimensional polytope PnHj with Hi as in Lemma El Now we 
can iterate this argument. If during this iteration the outer 1-radius of a polytope P' has 
to be computed, then by Proposition El the minimal enclosing 1-cylinder touches at least 
dim(P') + 1 affinely independent vertices. From the same iterative argument it follows 
that Rj(P) = R j+u - n -i(P n P) for some (u - l)-flat P. 

Suppose S = conv{fW, . . . , t>( n+1 )} is a simplex in E n , and J an [n — j)-flat, such that 

dist(w (1) , J) = ■■■ = dist(f ( "~- J+2) , J) = P^conv^ 1 ),...,^-^}) 
> dist(f ( ™"- ?+3) , J) > > dist(f (ri+1) , J). 

Then obviously Rj(S) = Pi (conv-fi/ 1 ) , . . . , t)( n ^'+ 2 )}) and at most n — j + 2 vertices are 
situated on the boundary of the minimal enclosing j-cylinder. 



The last point which remains to show is that every minimal enclosing j-cylinder of the 
regular simplex T n also circumscribes T n . Because of Proposition El it suffices to show 
that the value p* n in ()3.4j) is positive for all 1 < j < n — 1, showing that b) in Lemma El 
never holds for T n . 

However, in almost all cases the desired circumscribing property follows already from 
[Td] . see also |2j. 
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Proposition 7. For 1 < j < n it holds Rj(T n ) > y-^zf- Ifn is odd or j £ {l,n — 1}, 



then Rj(T n ) = d and every minimal enclosing j- cylinder of T n is a circumscribing 
j -cylinder of T n . 

We can easily apply Proposition [7| to compute p* n if n is even. 

Lemma 8. Let 1 < j < n — 1. If S = T n then always case a) in Theorem^ holds. 

Proof. We can assume j > 2, since otherwise b) cannot hold. We use the notation as in 
Lemma El Because of Proposition [7| it suffices to consider the case where n is even, and as 
mentioned above the proof is complete if we show that p* n is positive. Since Proposition [7| 
yields i? n _i(T n_1 ) = y/(n-l)/n, we have v { n +1) = ^2 - (n - l)/n = ^(n + l)/n. Also 
by Proposition [7[ i?j„i(T n_1 ) = — l)/n and therefore 

p* = n ~ J + 2 > o. 
2y/n(n + l) 

□ 

Choosing an optimal (n — l)-cylinder among those parallel to a facet of T n with p* n = 
gives an upper bound for the outer (n — l)-radius of a regular simplex, 

Rn-l(T U ) < 



2 x /n(n+l) 

2n - 1 



2y/n(n + l) ' 

Theorem 121 states that for even n this bound is tight. 

4. Reduction to an algebraic optimization problem 

In this section, we provide an algebraic formulation for a minimal circumscribing j- 
cylinder J + p(M n+1 D 7^q) of the regular simplex T n in standard embedding. Let J = 
p+linjs'- 1 ), . . . , s^ n ~^} with pairwise orthogonal s^ 1 ', . . . , s^ n ~^ G 5 n_1 , and p be contained 
in the orthogonal complement of lin{s < - 1 ' ) , . . . , s^ n ~^}. The orthogonal projection P of a 
vector z G H.™ onto the orthogonal complement of linjsW, . . . , s^"- - ^} (relative to Ti™) can 
be written as P(z) = (I — $^fc=i s^ k \s ( - k ^) T )z, where I denotes the identity matrix. Hence, 
for a general polytope with vertices . . . ,i/ m ) (embedded in H,™) the computation of 
the square of Rj can be expressed by the following optimization problem. Here, we use 
the convention 



min p 2 



(i) s.t. (p — Pv^) 2 < p 2 , i = l,...,m, 

(4.1) (ii) p- S W = 0, k = l,...,n-j, 

(iii) sW, . . . , s( n ~fi G iS™ -1 , pairwise orthogonal, 

(iv) p G HI- 
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In the case of T™, (i) can be replaced by 



n-j 



(4.2) (V) U_ e (i)+^ a (*) a (*)J = p \ i = l,..., n + l, 

where the equality sign follows from Theorem ^ By (ii) and G 5 n_1 , (i') can be 
simplified to 

n-j 

0") P 2 ~P 2 = E(^ ) ) 2 + 2 ^- 1 ' i = l,---,n + l. 

fc=i 

Summing over all i gives (n + l)(jo 2 — p 2 ) = (n — j) + 2 — (n + 1), i.e., p 2 — p 2 = We 
substitute this value into (i") and obtain p% = \ — Efc=i( s i^) 2 J • Hence, all the 

can be replaced in terms of the s\ , 

, ^ 2 (2 + (n-j))(2- (n-j)) l^f^ /( k),X j~l 

4(n+i) 4 trvfc / n+1 

1 n+1 n-j 

i=i fc'=i 

We arrive at the following characterization of the minimal enclosing j-cylinders: 

Theorem 9. Let 1 < j < n. A set of vectors s^ 1 ', . . . , s^ n ~^ G 5 n_1 spans the underlying 
(n — j) -dimensional subspace of a minimal enclosing j-cylinder ofT n C if and only if 
it is an optimal solution of the problem 

n+l /n-j N 1 

mm E ( E(«i ) 2 
i=i \fc=i 

s-t. EE(«i )*i = °> k = l,...,n-j, 

1=1 fc'=l 

sW,...,s' n_ ^ G 5 n_1 pairwise orthogonal. 

It is easy to see that in case j = n — 1 the program (|4.4|) reduces to (jl.lj) stated in the 
introduction. 

By (|4.3|h in order to prove i? n _i(T n ) = (2n — 1)/ (2yn(n + 1)) for even n, we have to 
show that the optimal value of (jl.lj) is 1/n. We apply the following statement from [3 . 

Proposition 10. Let n > 2. T/ie direction vector (si, . . . , s n+ i) T o/ an?/ extreme circum- 
scribing (n — l)-cylinder ofT n satisfies \{s\, . . . , s ra +i}| < 3. 

For completeness we repeat the short proof. 

Proof. We can assume n > 3. Let s G 5 n ~ 1 be the axis direction of a locally extreme 
circumscribing (n — l)-cylinder. Let f(s) := El=i s i De ^ ne objective function from (jl.lj) . 
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let g 1 ( 



n+1 

i 



1 «i, M^J 



s ? — 1> anc ^ S^ 5 ) := X^] 1 s «- necessary condition 



n+l 



for a local extremum is that for any pairwise different indices a,b,c,d G {1, . . . , n + 1}, 



det 



/ _ §L 
' ds a 

-2L 

9s b 

_«£ 

\ 

\ 9s d 



3s b 

dgi 
ds c 
dgi 
ds d 



892 
ds a 
9g2 
9s b 
9g2 
ds c 
9g2 
9s d 



9g3 
9s a 
9gs 
9s b 
893 
9s c 

9gz I 

9s d / 



-24 det 



V 



Sb 
S c 
Sd 



\ 



J 



The latter is a Vandermonde determinant, which implies \{s a , s&, s c , Sd}\ < 3. 



□ 



Using Proposition El can be written as the following polynomial optimization 

problem in six variables with additional integer conditions. 

min kis\ + k 2 s\ + k 3 s\ 

(i) s.t. kisl + k 2 s\ + k 3 s\ 

(ii) k v s\ + k 2 s\ + k 3 s\ 
(hi) fciSi + k 2 s 2 + k 3 s 3 
(iv) k\ + k 2 + k 3 

si,s 2 , s 3 el, k x ,k 2 ,k 3 



(4.5) 



0. 
1. 

0, 

n+l 



Since the odd case of Theorem |21 is well-known fUJ we assume from now on that 
n is even. The mindful reader will notice that for odd n the optimal value of (14.51) 
coincides with the optimal value of the real relaxation (where the condition ki, k 2 , k 3 G N 
is replaced by ki,k 2 ,k 3 > 0). 

For k 3 = the equality constraints in (|4.5j) immediately yield k\ = k 2 = (n+ l)/2 G" N, 
and similarly, for s 2 = s 3 we obtain k\ = k 2 + k 3 = (n + l)/2 N. Hence, we can assume 
that si, s 2 , and S3 are distinct and ki,k 2 ,k 3 > 1. Moreover, for S3 = the resulting 
optimal value is 1/n which will turn out to be the optimal solution. Finally, by (iii), not 
all of the Si have the same sign. Hence it suffices to show that for si < and s 3 > s 2 > 
every admissible solution to the constraints of (J4.5|) has value at least 1/n. 

The linear system of equations in k%, k 2 , k 3 defined by (i), (ii), and (iii) is regular, which 
is easily seen from the Vandermonde computation 



det 



4 

si 



S' 2 
6 2 



SlS 2 S 3 (si - s 2 )(si - s 3 )(s 2 - S3) ^ . 



s l S 2 S3 



S2 + S 3 



Solving for ki,k 2 , k 3 yields 

(4.6) ki 

(4.7) k 2 
(4-8) h 

S3[S 3 - Si)(s 3 - S 2 ) 

Since all factors in the denominators are strictly positive, (14. 7|) and (J4.8j) imply in partic- 
ular si + S3 > and si + s 2 < 0. 



-sx(s 2 - si)(s 3 - Si) 

Sl + gg 
s 2 (s 2 - si)(s 3 - s 2 ) ' 
-(si + s 2 ) 
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Using (iv) in ()4.5j) we can express one of the Sj by the others. Solving for s 2 gives 

s\ + s 3 



(4.9) s 2 



(n + l)sis 3 + 1 



Note that the denominator in (J4.9j) is strictly negative. 

Our main goal now is to show k± < (n + l)/2 and then to use the integer condition to 
deduce k\ < n/2. In order to achieve this, substitute ()4.9|) into the inequality Si + s 2 < 
which (in connection with si + s 3 > 0) allows to conclude S3 > sf > l/(n + 1). Then 
substitute ()4.9|) into ()4.6|) which yields 

n+ l (( n + i) s 2 _ + 1)s3 ( S3 _ Si) _ 2 ) 

M ^ — " 



2(s 3 - s 1 ){{n + l)s 2 l s 3 + 2s 1 + s 3 ) 
((n + 1)4 - l)(((n + 1)4 - 1) - ((n + l>i£3 + 1)) 
2(s 3 - si)(((n + l)sis 3 + l)si + (si + s 3 )) 



< 



since all factors within the last fraction are positive. Hence, k\ < (n + l)/2 and since it 
is an integer k\ < n/2. Similarly, although we do not need this, one can show k 3 < n/2. 
However, note that this bound does not hold for fc 2 . 
Now it follows from k\ < n/2 and (J4.6)) that 

< 2 — ns\ — 2s\ — n 2 s\s 3 + n 2 s\s\ — ns\s 3 + ns\s\ — 2ns\ + ns\s 3 
= -2{ns\ - l)((n + l)sis 3 + 1) + s 3 (si + s 3 )(n(n + l)s? -n-2). 

That means at least one of the two terms of the sum must be non-negative which gives 
s\ < —1/y/n. Moreover, s\ + s 3 > also implies s 3 > 1/y/n. 

Finally, we show that for any admissible solution to the constraints of (J4.5)) the objective 
value is at least 1/n. Replacing ki, k 2 , k 3 and s 2 in the objective function via (|4.fij) - (|4.9j) 
and using the inequalities — s\ < s 3 and s\ < —1/y/n obtained above we get 

k lS t + k 2 si + k 3 st = + + 

1 2 3 n + l ( n + !)(_(„ + i) SlS3 _i) 

> 1 | ((n + 1)4-1) > 1 | 1 = 1 
— n + l n + l — n + l n{n + 1) n 

Hence, the optimal value of (|4.5jl is 1/n, and by our remark before Proposition HTH this 
completes the proof of Theorem 



5. Connections to the Positivstellensatz 

We close the paper by discussing the greater difficulty of the even case of computing 
i? n _i(T n ) compared to the odd case, by analyzing problem (jl.lj) from the viewpoint of the 
Positivstellensatz j20j. This theorem states the existence of a certificate whenever a system 
of polynomial equalities and inequalities does not have a solution, and it can be regarded 
as a common generalization of Hilbert's Nullstellensatz and of linear programming duality. 
For our purposes, it suffices to consider the following version of Putinar (see [El EH])- For 
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n G N let M[x] = M[xi, . . . , x n ] denote the ring of polynomials in xi, . . . , x n , and let 

^M[x] 2 = 6f for some Jfe e N, &i, . . . A e R[x) J 

be the set of all finite sums of squares of polynomials. Set g := 1. 
Proposition 11. If the polynomials f,gi,...,g m G ~R[x] satisfy 

f(x)>0 for all x G S := {x G R n : g x {x) > 0, . . . , g m (x) > 0} 

and 

(5.1) M : = : a ,...,a m G ^M[x] 2 j 

contains the polynomial 1 — EILi x h then f G M. 

Hence, in order to prove that the optimal value of (jl.l|) is bounded from below by some 
given value a, it suffices (by compactness of the admissible set) to show the existence 
of such a representation for f(s) := YH=i s t ~ a + £ 1X1 terms of gi(s) := ^2^=1 s h 
92(s) := -£r=i^ 9s(s) := EIL + M " 1, := -ES*? + 1, ffsW := ES 

g 6 (s) := - Y,7=i s i for ever y £ > °- 

Bounding the degrees of the polynomials cx^ by a fixed constant in (|5.1|) serves to give 
lower bounds on the minimum. These relaxations can be computed by semidefinite pro- 
gramming (SDP) and are at the heart of current developments in SDP-based constrained 
polynomial optimization (see HHj). 

For the case n odd of (jl.lj) there exists a simple polynomial identity 

n+l 1 „ / n+1 \ n+1 , 1 \ 2 

i=i \i=i / i=i v 7 

which shows that the minimum is bounded from below by l/(n + 1), and since this 
value can be attained by s\ — . . . — S(rt+i)/2 = _ S( n +3)/2 = ■ ■ ■ = — s n +i = 1/y/rT+T, 
the minimum is l/(n + 1). For any e > 0, adding e on both sides of (|5.2jl yields a 
representation of the positive polynomial on the left side as a sum of squares of the <?j. 
Note that for every odd n this representation uses only polynomials aigi of (total) degree 
at most 4. 

For the case n even (with minimum 1/n) the situation is quite different. A computer 
calculation using the Software GloptiPoly jTU] showed that already for n = 4 it is 
necessary to go up to degree 8 to find the Positivstellensatz-type certificate of optimality. 
Since from a practical point of view the computational efforts of this calculation drastically 
increase with the number of variables, we do not know up to which degree it is necessary 
to go for n — 6. 

Acknowledgments. We would like to thank Keith Ball for pointing out a substantial 
simplification in the proof of Theorem El and Pablo Parrilo for helpful discussions. 
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